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1. Introduction 

This paper is the third in a series of four reports on our calculation of Bk using HYP-smeared 
staggered fermions. In the previous two reports we presented the results of our analysis using, 
respectively, SU(3) and SU(2) staggered chiral perturbation theory (SChPT) [HI 01. Here we explain 
how we estimate the finite-volume error, and discuss issues concerning the statistical error. The 
calculations use the MILC ensembles (C3) and (C3-2), the parameters of which are given in Table I. 



parameter 


value 


P 


6.76 (Nf = 2+1 unquenched QCD) 


l/a 


1.588(19) GeV 


geometry 


20 3 x 64 (C3) and 28 3 x 64 (C3-2) 


configurations 


671 (20 3 ) and 274 (28 3 ) 


measurements per conf. 


9 (20 3 ) and 8 (28 3 ) 


sea quark (asqtad) masses 


ami = 0.01, am s = 0.05 


valence quark (HYP-smeared) masses 


0.005, 0.01, 0.015, . . . , 0.04, 0.045, 0.05 



Table 1: Ensembles used to study finite volume effects. 



2. Finite Volume Correction: Theory 

Finite volume (FV) effects are predicted by ChPT, and are caused by pseudo-Goldstone bosons 
(PGBs) in loops propagating from their starting point to a periodic image. Experience indicates that 
the one-loop calculation is a good guide to the sign and rough magnitude of the FV effect (when 
compared to numerical results or higher order calculations), although it tends to underestimate 
the size by factors of two or so. We have calculated the one-loop FV effects using SChPT for 
both SU(2) and SU(3) formulae. These enter in the standard way through the chiral logarithmic 
functions: 



(2.1) 



£(X) =X[log(X/^ R ) + Sr(X,L 

1{X) = = - lo ^ X /^R) - 1 + 5 3 FV (*^) (2-2) 

where X = M 2 is a PGB mass-squared (in physical units), and the finite-volume corrections are 
given in terms of modified Bessel functions by 

8r(X,L) = 4 -Y^± (2.3) 

5 3 FV (X,L) = 2£tf (|H|z). (2.4) 

Here z = ML, with L the spatial box size in physical units, and n is a vector labeling the spatial 
image positions, taking values (1,0,0), (0, 1,0) etc.. 1 

1 Strictly speaking, we should also include images in the time direction as well, but, since L, S> L for the MILC 
lattices, these contributions turn out to be negligible compared to the errors in Bk- 
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When evaluating the image sums we find that keeping up to n 2 = 12 in 8™ and n 2 = 16 in 8™ 
gives sufficient accuracy. As a side note, for our range of z, 2.5 < z < 10, it is sometimes useful to 
have an approximate form to speed up the calculation. We find that 

5f v (z) » ^V 3/2 (26.0095 + 37.8217/z+ 112.694/z 2 -59.9524/z 3 ), (2.5) 
S 3 FV (z) « e~V 1/2 (12.8717 + 10.7318/z + 88.9023/z 2 -20.8847 /z 3 ). (2.6) 

gives an accuracy of better than 1 part in 10 3 over this range. Note that these are not asymptotic 
expansions — these forms fail for both small and large z. 

In our fitting to SChPT forms (as described in the companion reports ^ ||]) we have, so far, 
used the infinite-volume expressions. Thus it is important to check that FV corrections are small. 
Using the coefficients from the fit functions, we can calculate the one-loop expectation for the FV 
corrections for each fit. For the SU(3) fits, the largest contribution to the FV corrections comes 
from terms induced by discretization and truncation errors. These are not well determined, with 
the result that different fits differ even in the sign of the predicted FV effect. Thus we can at most 
use these calculations to estimate the magnitude of the FV effect. 



FIT N-T2: FV shift in BK 




-0.004 - 



Figure 1: Finite volume shift in B K as predicted by SU(3) SChPT for the 20 3 x 64 coarse MILC lattice, using 
the "N-T2" fit [^J. The horizontal axis orders the 55 valence "kaons" into groups in which am y is fixed and 
am x is increased from 0.005 up to am y . The first group (10 leftmost points) have am y — 0.05 « amf hys , the 
next group (9 points) have am y — 0.045, etc.. The rightmost point has am x = am y = 0.005. The extrapolation 
to the physical kaon mass is dominated by the leftmost group. Note that the lowest "row" of points all have 
am x = 0.005. 

As an example, we show in Fig. p] the FV correction for the ensemble (C3) using the fit-type 
which gives rise to the largest FV correction. We see that the effect depends dominantly on the 
light-quark valence mass m x , and much more weakly on the strange valence mass m y . It grows 
rapidly as m x decreases, which is expected given the exponential dependence on z. The size of the 
FV effect is, however, small compared to Bk ~ 0.5 — the relative contribution does not reach 1%. 
We also note that on the larger, 28 3 , lattices, the FV correction is smaller by almost an order of 
magnitude, and thus negligible compared to other errors. 

The situation with SU(2) SChPT is shown in Fig. ^, here with results on both coarse and 
fine lattices. The NLO prediction is more reliable in this case, since the coefficient of the chiral 
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FV effects in BK with SU(2) ChPT on ensemble (C3) 



FV effects in BK with SU(2) ChPT on ensemble (Fl) 
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Figure 2: Finite volume correction to Bk from SU(2) SChPT as a function of the squared mass of the pion 
composed of light valence quarks, M\ =X in GeV 2 . Left panel: coarse ensemble (C3); right panel: fine 
ensemble (Fl). The lower limit of the curves corresponds to the minimum values of M\ in the simulations. 

logarithms is well determined by the fits. We see again that the FV effect grows rapidly as m x <=< M\ 
decreases. The size of the effect is, however, very small, significantly smaller than predicted by 
the SU(3) fit. One peculiar feature is that the sign changes (and magnitude increases) moving from 
coarse to fine lattices (which are otherwise similar). This is a reflection of the fact that, as one 
approaches the continuum limit, taste splittings reduce, more pions become lighter, and FV effects 
increase. In this specific case there are various canceling contributions and it turns out also that the 
sign flips. 

We draw two main conclusions from these results. First, the FV correction to Bk is small, 
comparable to or smaller than the statistical error for our most physical "kaon" (an error which is 
0.6% on the coarse ensemble). Thus we are unlikely to make an error of larger than « 1% by 
leaving FV corrections out of the fits. Second, we do not know in detail the expected form, or even 
the sign, of the FV correction. Note that the true FV correction for given valence masses has a 
definite value, one that SU(2) and SU(3) ChPT should agree upon. The disagreement between the 
two estimates is thus an indication that the NLO predictions are quantitatively unreliable. 

This topic deserves further study. In particular we intend to include one-loop FV corrections 
in future fits. In the meantime we turn to a more direct method for estimated the FV error. 

3. Finite Volume Correction: Direct Measurement 

In order to study the finite volume effect numerically, we compare results on the two coarse 
ensembles whose parameters are listed in Table |]. The lightest valence pion on the smaller lattice 
has m n L m 2.68. This is below the rule-of-thumb minimum "safe" value of 3-4, but we expect 
to have more leeway when calculating kaon properties (such as Bk), since the pion only enters 
through loops. Indeed, the estimates of the previous section suggest that FV effects remain small. 
We also note that the same estimates indicate that FV effects are almost negligible on the larger, 
28 3 , lattices. 

We have adjusted the number of measurements per configuration so that the statistical weights 
of the two ensembles are very similar: 



ratio of stat. weights 



20 3 x671 x9 



= 1.004. 



(3.1) 



28 3 x 274 x 8 
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Figure 3: One-loop matched Bk vs. m\ for the 10 degenerate combinations (left) and for all 55 degenerate 
and non-degenerate combinations (right). 



In Fig. 51 we show the results for Bk on the two lattices, and find that the errors are indeed com- 
parable. The results on the larger lattice lie systematically higher, although for each point the 
difference is only about 1 a. We fit both data-sets using our preferred Bayesian "N-BT7" fit for the 
SU(3) analysis, and the "4X3Y-NNLO" fit for the SU(2) analysis. 2 The SU(2) "X-fits" are shown 
in Fig. | 
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Figure 4: One-loop matched Bk vs. Xp for the SU(2) analysis on the 20 3 volume (left) and on the 28 3 
volume (right). Here Xp is the mass-squared of the Goldstone taste "pion" composed of light valence quarks 
of mass am x . See Refs. [|lj ^| for more details. 

The results of the fits are summarized in Table ||. The difference between Bk on the two vol- 

I- 



2 The details of these fits are discussed in Refs. 
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umes is 2% for the SU(3) analysis, which is about twice the statistical error, and thus probably 
significant. For the SU(2) analysis the finite volume effect is smaller, 0.9%, and, being compa- 
rable to the statistical error, less significant. We interpret the difference between SU(3) and SU(2) 
analyses as being due to greater simplicity of the SU(2) fit form. The size of these effects are 
not unexpected given the theoretical results of the previous section. We use these differences as 
estimates of the FV systematic error in the respective fits, as quoted in Refs. ^]. 



volume 


SU(3) 


SU(2) 


20 3 


0.5599(57) 


0.5645(49) 


28 3 


0.5730(57) 


0.5694(48) 



Table 2: Results for Bk using SU(3) and SU(2) staggered ChPT. The values quoted are in the NDR scheme 
at a scale of fi = 2 GeV. 



4. Increasing the Statistics and Autocorrelations 

It is well know that one can obtain more information from each gauge configuration by using 
multiple time-positions for the sources. An important issue, however, is the extent to which these 
measurements are independent. We have tested this by using 9 randomly chosen "starting times" 
for the wall sources on each configuration in the coarse ensemble (C3), each with a different random 
seed for the noise in the source. We average the raw data from these 9 sources, and then compute the 
error using single elimination jackknife. We compare the results so obtained (labeled "x9") to those 
obtained using a single measurement/configuration ("xl") in Fig. ||, displaying only degenerate 
combinations for clarity. We observe that the errors are reduced by a factor of 3, indicating that the 
different measurements on each configuration are almost independent. 
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Figure 5: Bk vs. m\ on ensemble (C3) for one measurement/conf. and for 9 measurements/conf. 
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We have also studied autocorrelations between configurations by looking at the dependence of 
the statistical errors on the bin size (i.e. the number of configurations we bin together) In Fig. ^, 
we display the result for the error on the 2-point correlator 



C(t) = (A 4 (t)P(0)) 



(4.1) 



(where both operators have the Goldstone taste ^5) at t = 10. For the xl data, the error in the 
error bar is too large to determine whether there is an autocorrelation. 3 For the x9 data, however, 
we do observe about a 20% increase between a bin size of 1 and 3-6. Fortunately this is a small 
effect. Nevertheless, it is clear that having higher statistics allows us to better understand the 
autocorrelations, and we intend to increase the statistics on other ensembles as well. 
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Figure 6: Error in C(t = 10) [see text] vs. bin size for the xl data (left) and the x9 data (right). 
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Details of how we estimate the error in the error bar will be given in Ref. 



